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Решена пространственная нелинейная нестационарная задача теплопроводности для составного параллелепипеда с локально 
сосредоточенными источниками тепловыделения и неоднородными граничными условиями. Использован метод конечных раз- 
ностей и неравномерная разностная сетка. Сделан вывод о существенной неоднородности температурных полей в сечениях, 
проходящих через источники тепловыделения. 


Одной из актуальных проблем настоящего вре- 
мени является проблема теплоэнергосбережения, 
решение которой невозможно путём эмпирическо- 
го анализа различных технологических, проектных, 
материаловедческих и конструкторских решений и 
предложений. Эффективным инструментом поиска 
решений проблемы теплоэнергосбережения может 
быть только математическое моделирование комп- 
лекса процессов, протекающих в реальных систе- 
мах, потребляющих тепловую энергию, с последую- 


щей опытной проработкой наиболее привлекатель- 
ных решений и схем. Но до настоящего времени не 
опубликованы результаты исследований по созда- 
нию теоретических основ процессов простран- 
ственного нестационарного теплопереноса в систе- 
мах-потребителях тепловой энергии. Целью данной 
работы является решение задачи пространственно- 
го нестационарного теплопереноса в объекте, 
представляющим собой замкнутый объем с локаль- 
но сосредоточенными источниками тепловыделе- 
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ния и неоднородными граничными условиями на 
внешних и внутренних границах области решения. 



Рис. 1. Область решения рассматриваемой задачи: 1) кир- 
пичная стена ; 2) воздух ; 3) стеклянное окно ; 4) дере- 
вянный пол; 5) бетонная плита; 6) локально сосредо- 
точенные источники тепловыделения 


Рассматривается краевая задача теплопровод- 
ности для области, представленной на рис. 1. Об- 
ласть решения включает 12 параллелепипедов, 
имеющих разные размеры и теплофизические ха- 
рактеристики. На границах между всеми паралле- 
лепипедами и на границах с внешней по отноше- 
нию к рассматриваемому объекту средой выставля- 
лись соответствующие граничные условия. Иссле- 
дуемая в данной работе модель объекта теплопот- 
ребления (рис. 1) достаточно хорошо описывает 
помещение, отапливаемое за счёт централизован- 
ной или местной системы отопления в условиях 
интенсивного теплообмена с внешней средой. 

Задача решена в рамках модели теплопровод- 
ности без учёта возможного теплопереноса за счёт 
естественной или вынужденной конвекции. В та- 
кой постановке процесс переноса тепла в анализи- 
руемом объекте описывается системой нестацио- 
нарных уравнений теплопроводности с нелиней- 
ными граничными условиями. 

Математическая модель. Распределение тепла в 
рассматриваемой области описывается уравнением 
теплопроводности [1]. 
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где Т=Т(і,х,у,і); Т ; - температура; х,у,і- декартовы 
координаты; с, - удельная теплоемкость; р, - плот- 
ность; А, - теплопроводность /-ого материала; / = 1 - 
кирпичная стена; 2 - воздух; 3 - стеклянное окно; 
4 - деревянный пол; 5 - бетонная плита; Г - время. 

Начальные и граничные условия для сформули- 
рованной задачи имеют вид: 

Начальное условие: 

Т((,х,у,г)[ =0 =Т 0 =сошД. 


Граничные условия: 

- на грани, разделяющей внешнюю среду и рас- 
сматриваемый объем 

-А. — = а(Г - Т.) + ео(Т* -Т. 4 ), 

і -ѵ х е і ' 4 е і 

оу 


где і = 1,3,5; поскольку на данной грани находятся 
такие материалы как: / = 1 - кирпич; 3 - стекло; 
5 - бетон. 

- на всех участках области решение, где происхо- 
дит сопряжение материалов с различными теп- 
лофизическими параметрами, заданы условия 
4-го рода 

Т = Т. і = 1 , 5 

ВДе 

' Эх* 2 дх к к = 13 


- на всех внешних границах рассматриваемого 
объема, кроме грани, на которой осуществляет- 
ся теплообмен с внешней средой, заданы усло- 
вия симметрии 



= 0 , где 


і=1,5 

к=Ъз’ 


- на локально сосредоточенных источниках теп- 
ловыделения заданы условия 1-го рода 

Т =Т., = СОП8І, 

Ьаі ’ 


здесь а - коэффициент теплообмена между внеш- 
ней средой и рассматриваемой областью решения; 
Т е - температура внешней среды; е - приведенная 
степень черноты; Г 0 - начальная температура рас- 
сматриваемого объекта; Т ш - температура на ис- 
точниках тепловыделения; а - постоянная Стефа- 
на-Больцмана. 

Предполагалось, что источники тепловыделе- 
ния имеют постоянную в течение всего времени 
температуру, а на их границах выполняются гра- 
ничные условия первого рода. 

Сформулированная краевая задача решена ме- 
тодом конечных разностей [2]. Особенность рас- 
сматриваемой задачи состоит в том, что область ре- 
шения включает несколько элементов с сущест- 
венно различными по величине теплофизически- 
ми характеристиками. При этом самая большая по- 
добласть решения имеет размеры в сотни раз отли- 
чающиеся от размеров характерных подобластей с 
существенно отличными коэффициентами теп- 
лопроводности и температуропроводности. В связи 
с этой особенностью решение сформулированной 
краевой задачи невозможно с применением регу- 
лярной и равномерной разностной сетки в задан- 
ном диапазоне изменения основных параметров. 
Поэтому применялась неравномерная разностная 
сетка со сгущением вблизи границ раздела подоб- 
ластей с различными теплофизическими характе- 
ристиками материалов. Для решения разностных 
аналогов исходных дифференциальных уравнений 
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использовался локально-одномерный метод, а для 
решения одномерных разностных уравнений - 
итерационная неявная разностная схема [2, 3]. Вы- 
бор значений сеточных параметров проводился на 
множестве сеток исходя из условий обеспечения 
устойчивого и сходящегося решения при мини- 
мально допустимом числе узлов разностной сетки. 



ИСТОЧНИКИ 

те п л овы д е л е и ия 


Рис. 3. Сечение области решения плоскостью г = 0,7 м 


Рис. 2. Распределение темпера туры в сечении г = 0,7 м в мо- 
мент времени 1 = 32400 с. Точки по осям х и у — ус- 
ловные координаты, отвечающие номерам узлов 
разностной сетки области решения 

На рис. 2 приведены типичные результаты ре- 
шения сформулированной краевой задачи в виде 
температурного поля Т(х, у ) в фиксированном се- 
чении і = 0,7 м в момент времени { = 9 ч. Это сече- 
ние соответствует по конфигурации границ рис. 3. 


Численный анализ проведён при следующих 
значениях: а = 10 Вт/(м 2 -К), е = 0,3, Т е = 258 К, 
Т 0 = 293 К, Т ш = 333 К. Длина комнаты -5 м; ши- 
рина -4 м; высота - 3 м; толщины стен - 0,3 м. 
Длина окна -2 м; высота -1,5 м; толщина стекла - 
0,005 м; расстояние между стеклами -0,1 м. Длина 
батареи - 0,5 м; толщина - 0,1 м; высота - 0,6 м. 
Толщина пола - 0,2 м; расстояние от пола до бата- 
реи - 0,2 м; расстояние от батареи до окна - 0,2 м; 
расстояние от окна до потолка -0,5 м; толщина по- 
толка - 0,2 м. Расстояние от боковой стены до ок- 
на-1,5 м; расстояние от боковой стены до первой 
батареи и от второй батареи до стены -0,5 м. 

На рис. 2 хорошо видна существенная неодно- 
родность температурного поля, характерная для се- 
чений, проходящих через источники тепловыделе- 
ния. В других сечениях (рис. 4, а) эта неоднород- 
ность проявляется менее отчетливо, но сохраняет- 
ся и на достаточно большом удалении от источни- 
ков (рис. 4, б). 

Следует особо отметить, что для решения одно- 
го варианта рассматриваемой краевой задачи тре- 
буются значительные ресурсы ПЭВМ. Исходя из 
того, что известны решения одномерных нестацио- 





Рис. 4. Распределение температуры в момент времени I = 32400 с в сечении: а) г = 1,1 м; б) г = 1,95 м 
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нарных краевых задач теплопроводности для ана- 
логичных объектов исследования [4, 5], проведены 
сравнения результатов вычислений полей темпера- 
тур в идентичных условиях для одномерной и 
пространственной постановок. На рис. 5 приведе- 
ны результаты сравнения для сечений і= 1,95 м, 
х = 2,8 м. 



Рис. 5 . Распределение температуры при одномерном и прост- 
ранственном случаях в момент времени I = 32400 с 

Хорошо видно, что отличия по значениям Т сос- 
тавляют 1,5. ..2, 0 К. Сравнение среднеобъёмных 
значений температуры, вычисленных при исполь- 
зовании одномерной и пространственной моделей, 
приводит к аналогичным результатам. Учитывая 
специфику проблемы теплоэнергосбережения, 
можно сделать вывод [6], что решение простран- 
ственной краевой задачи приводит к существенно 
более точным оценкам значений температур внутри 
типичных объектов теплоснабжения. Следует отме- 

СПИСОК ЛИТЕРАТУРЫ 

1. Лыков А.В. Теория теплопроводности. — М.: Высшая школа, 
1967.-600 с. 

2. Самарский А.А. Теория разностных схем. — М.: Наука, 1977. — 
656 с. 

3. Березин И. С., Жидков Н.П. Методы вычислений. — М.: Физ- 
матгиз, 1962. — Т. 2. —640 с. 

4. Иванов В. В., Карасева Л.В., Тихомиров С.А. Моделирование 
процессов теплопереноса в многослойных ограждающих 
конструкциях // Труды третьей Российской национальной 
конференции по теплообмену. — М.: Издательство МЭИ, 2002. 
-Т. 7. -С. 131-134. 


тить, что рассматриваемая постановка может быть 
существенно усовершенствована путём учёта меха- 
низма теплообмена внутри рассматриваемого объ- 
екта за счёт естественной конвекции. В этом случае 
результаты анализа температурных полей могут 
быть уточнены. При реализации же одномерных 
постановок для такого уточнения необходима экс- 
периментальная информация по коэффициентам 
теплообмена на всех внутренних границах области 
решения. Определение значений этих коэффици- 
ентов представляет достаточно сложную задачу, для 
решения которой необходимы серьёзные усилия в 
каждом конкретном случае. В то же время, при реа- 
лизации пространственной модели необходимость 
проведения экспериментов для определения исход- 
ных данных отпадает. Процесс моделирования мо- 
жет быть проведен в этом случае и для сложных 
конфигураций исследуемых объектов при наличии 
сведений о теплофизических характеристиках мате- 
риалов, используемых в данном объекте. Эти харак- 
теристики хорошо известны не только для традици- 
онных строительных материалов (кирпич, бетон, 
железобетон, древесина, стекло) [7], но для новых 
конструкционных материалов (пластмассы, тепло- 
изоляционные материалы, пенобетоны и др.). 

Решена пространственная нестационарная за- 
дача теплопроводности для замкнутого объема с 
локально сосредоточенными источниками тепло- 
выделения, проведены сравнения температурных 
полей для пространственной и одномерной поста- 
новок рассматриваемой задачи. Представленные 
результаты позволяют делать заключение о неод- 
нородном распределении температуры. 
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